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Assuming that the effect of the residuai interaction beyond mean-fieid is weak and has a short 
memory time, two approximate treatments of correlation in fermionic systems by means of Marko- 
vian quantum jump are presented. A simplified scenario for the introduction of fluctuations beyond 
mean-field is first presented. In this theory, part of the quantum correlations between the resid- 
ual interaction and the one-body density matrix are neglected and jumps occur between many-body 
densities formed of pairs of states D — |<3? a ) / ($(, | $ a ) where |<3> a ) and |$(,} are antisymmetrized 
products of single-particle states. The underlying Stochastic Mean-Field (SMF) theory is discussed 
, and applied to the monopole vibration of a spherical 40 Ca nucleus under the influence of a statistical 

■ ensemble of two-body contact interaction. This framework is however too simplistic to account for 

both fluctuation and dissipation. In the second part of this work, an alternative quantum jump 
£\J , method is obtained without making the approximation on quantum correlations. Restricting to two 

particles-two holes residual interaction, the evolution of the one-body density matrix of a correlated 
system is transformed into a Lindblad equation. The associated dissipative dynamics can be simu- 
lated by quantum jumps between densities written as D — [<&) (<£>| where $) is a normalized Slater 
determinant. The associated stochastic Schroedinger equation for single-particle wave- functions is 
given. 
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The description of quantum self-interacting systems with many degrees of freedom is common to many fields of 
physics, including Bose-Einstein condensates, atomic clusters and nuclear systems. A strikingaspect related to this 
problem is the emergence of well ordered motion at the same time as complexity and chaos In many situations, 

the self-consistent mean- field theory provides a suitable framework to describe ordered motions. It also corresponds 
to one of the most useful method to study the static and dynamical properties of self- interacting systems [310, 0. 
' However, it often turns out that mean-field reproduces average properties of one-body observables but underestimates 
dissipative and fluctuating aspects. This can directly be assigned to the absence of two-body effects beyond mean- 
field. In the nuclear context, it has been proposed to extend mean-field theory by considering that one-body degrees 
rj-( of freedom represents a subsystem which is coupled to more complex internal degrees of freedom. By doing so, the 
• • ■ problem of self-interacting system has been mapped to an open quantum system problem. In nuclei, extensive work 
has been devoted to the formal derivation of dissipative quantum mechanics Rl and/or related stochastic equations 
for fermions, including Markovian and non-Markovian effects [7llci IsL IToL [Til \12L fl3L Il4j ]. These approaches have in 
common that the residual part of the interaction introduces disorder on top of the mean field. These theories end 
with rather complex transport equations which are hardly applicable in realistic situations 15|. In fact, only recently, 
the theory proposed in ref . [Tj| has been applied to small amplitude collective vibrations [lq . Its application to large 
amplitude motion in non-equilibrated quantum many-body dynamics remains an open issue |l5l Il7| . 

During the past decades, large theoretical efforts have been devoted to the development of Monte-Carlo methods to 
describe the static properties of many-body interacting systems |18| . Recent applications to nuclear physics have shown 
that stochastic methods can successfully be applied to describe the structure of nuclei 0] . These methods can also 
be extended to dynamical problems [BJ and has been used recently to treat the exact dynamics of interacting bosons 
pof or fermions [2l| in schematic cases. As underlined in ref. [2^, it is highly desirable to provide approximate theory 
in order to describe dissipation in many-body systems taking advantage of recent advances in Monte-Carlo methods. 
The present work is an exploratory study devoted to the description of dissipation and fluctuations in dynamical 
problems using quantum jump techniques. Starting from a perturbative treatment of the residual interaction, two 
strategies are used to transform the many-body dynamics into a stochastic process. In the first strategy, neglecting 
part of the quantum correlations, a stochastic theory that may be economical in terms of numerical implementations 
is obtained. While appropriate to treat fluctuations beyond mean-field, it is however not suitable for dissipation. 
In the second part of this article, we show that the approximate treatment of quantum correlations can be avoided 
leading to a more general framework. In that case, the dynamical equation of motion of one-body degrees of freedom 
can be mapped into a Lindblad equation |23l l24j generally found in the theory of open quantum systems plj | . The 
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associated jump process is finally discussed. 



A. Generalities on perturbation theory and stochastic mechanics 



We consider a many-body fermionic system described by a two-body Hamiltonian H. We assume that the mean- 
field theory already provides a good approximation of its static and dynamical properties. In this case, the N-body 
wave function can be replaced by an antisymmetrized product of single particle states interacting through an effective 
self-consistent mean- field, denoted by Hmf- Let us assume that the system is initially a Slater determinant, denoted 
by and let us introduce the mean-field propagator UMp(t' , t). The great advantage of mean-field theory, is 

that the propagated many-body state | = Umf^', t) \&(t)) remains a Slater determinant and the dynamical 
evolution of the system reduces to the evolution of its single-particle components. Thus, the many-body density D 
is approximated by D ~ |$(f')) Accordingly, all the information on the system is contained in the one-body 

density matrix, denoted by p, whose matrix elements are defined by (j \p\i) = TrfafajDj = (afdj). Mean-field 
theory does simplify the dynamical description of many-body systems by reducing significantly the number of degrees 
of freedom to follow in time. 

In nuclear physics, mean-field is often adequate to describe average properties of one-body observables but fails to 
account for fluctuations. At the wave-function level, this corresponds to a deviation of the mean-field trajectory from 
the exact dynamics \26L l27j . 

In order to illustrate this effect, we denote by 8v\ 2 the residual two-body interaction defined through 8v\% = 
H — Hmf- In the weak coupling regime, this deviation can be treated in perturbation theory and the state at time 
t' reads [l^ 



|*(f)) = Mt'))-^Jsv 12 ( S )m S ))ds~^T^J J 5v 12 ( S ')Sv 12 ( S )d S 'd S \<S>( S )) 



(1) 



Svi2(s) corresponds to the residual interaction written in the interaction picture, Svi 2 (s) = U MF (s, t)Svi2(t)UMF(s, t). 
While at the mean-field level, \^(t')) is replaced by |$(t')), due to the accumulated effect of Sv±2 in time, it is expected 
that the exact state becomes a more and more complex superposition of a large number of Slater determinants. 
Accordingly, the information on the system cannot be reduced anymore to the knowledge of the one-body density 
matrix only and higher order correlations must be accounted for. A natural extension of mean-field theory is to enlarge 
the number of degrees of freedom considered. This is done for instance in the Time-Dependent Density Matrix theory 
(TDDM) where both the one-body density matrix and the two-body correlation operator are followed in time |28| . 
These theories are however rarely applied due to the large number of degrees of freedom to consider |29j . 

An alternative way to account for correlations beyond mean-field is to use stochastic methods. Several strategies 
have bee n prop osed to introduce noise on top of the mean-field, either based on statistical ensemble of one-body 
densities [JEI El > random two-body interactions or phase shift [3, Q , or directly from Fermi golden rule 0, 0, Hif . 
The goal in all cases, is to simulate the dynamics of correlated system by averaging over an ensemble of mean-field 
trajectories. The great advantage in that case, is that the number of degrees of freedom to follow along each path is not 
increased compared to standard mean-field. However, these theories are rather complex and methods for numerical 
implementations are still missing. 

The aim of the present work is to discuss again the possibility to replace correlated dynamics by quantum jumps 
in the Hilbert space of Slater determinants. We use equation Q as a starting point and we restrict the discussion 
to the Markovian limit. In order to illustrate this hypothesis, we follow ref. We assume that the residual 

interaction induces random transitions treated as a statistical ensemble of two-body interactions acting on top of the 
mean-field. Eq. is then replaced by a set of evolutions with the same initial state and mean-field but with different 
residual interactions. We assume that the two-body operator has a gaussian distribution with a mean value Svi 2 = 
and a second moment denoted by 8v\ 2 . Here, the average is taken over different values of 8v\2- In nuclear systems, 
the residual interaction is expected to induce transitions on a shorter time scale (called correlation time and denoted 
by t) than the time associated to the mean field evolution (denoted by r re ;)0, l30|. r is related to the average 
autocorrelation function 8v\2{s')8vi2{s) which is approximated by 

8v 12 (s')Sv 12 (s) cx ^(s) e -(s-s') 2 /2^ (2) 

Using this approximation, we consider a time-scale At much larger than the time r but smaller than r re ;. In the 
following, this limit will be called "Markovian" or "short memory time" approximation. 

Using approximation , we consider in the following two limits for which the wave function evolution as given by 
can be replaced by quantum jumps between Slater determinants. The first case is a simplified scenario where part 
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of quantum correlations between 8vi 2 and p are neglected along the path. In this case, it is shown that the evolution 
can be formulated in terms of quantum jumps between many-body densities formed of pairs of Slater determinants. 
In a second part, we show that the perturbative dynamics can still be transformed into a quantum jump process even 
if quantal correlations are not neglected. In both cases, all equations necessary for applications are given explicitly. 



II. FLUCTUATIONS BEYOND MEAN-FIELD IN A SIMPLIFIED STOCHASTIC SCENARIO 

Let us first consider the perturbative evolution of an initial Slater determinant. Assuming equation J2J and the 
short memo ry tim e approximati on, the mean-field does not change over At. Then, the average evolution of the state, 
denoted by A |<£) = \V(t + At)) - |$(i)>, reduces to 

At„ „ rAt 



A |*) = —Bmf |*(t)) - w bv\ 2 |$(t)) . (3) 

This expression can also be regarded as an average over Markovian stochastic processes in many-body wave-functions 
space. To give a deeper insight we define the ensemble of antisymmetrized two-body residual interactions as 

£«i2(o0 = 4 X] a " a P ( a/3 \ Svi2 ( a )\ <*7) a i a s- (4) 

Here, a = N where all components are stochastic variables sampled according to gaussian probabilities with 

mean zero and ai 2 — 1. The number N of stochastic components defines the complexity of the process. The definition 
|3| includes the force proposed in ref. |9(. Eq. can be interpreted as the average over the quantum diffusion 

A |*) = {^-Hmf + ABSv 12 + \ (AB5v 12 fX |$(t)) , (5) 



where AB = iyrAt/Ti. In the following, we will consider last expression as a differential stochastic equation in Hilbert 
space |2a, |31| . We use the notation dB instead of AB and use the Ito rules of stochastic calculus [32j . Due to the 
two-body nature of 5vi 2 , eq. JSJ induces complex reorganization of single-particle degrees of freedom. After the jump, 
the state is not a priori a single Slater determinant. For applications, it is highly desirable to preserve the simple 
initial form of the state along the stochastic path. This could be achieved by invoking additional approximations 
described below. Following ref. [20L l2ll , we consider an initial density 

D = MJM (6) 

where |<&„) = A (Hi \cti)) and ($(,) = A (H; are two non-orthogonal Slater determinants formed of products of sin- 
gle particle wave packets denoted respectively by \cti) and |/3,-). The notation A(.) corresponds to the antisymmetrized 
product. We assume that both states follow the diffusion process described by eq. (JSJ) but with two independent 
sets of gaussian stochastic variables, denoted respectively by cr a and ov This case will be referred in the following 
as the " uncorrelated noise". The use of different sets of stochastic variables is at variance with standard quantum 
Monte-Carlo procedures that simulate density evolution given by Lindblad equations 25] . However, this assumption 
has been shown to be crucial in order to describe the exact dynamics of interacting systems with stochastic methods 

MM. 



A. Approximate stochastic mechanics in one-body density matrix space 

To approximate the diffusion process, we first focus on single-particle degrees of freedom. Under the approximation 

(afajSvlz) ~ (afaj) (5v 2 l2 ) + 2 (a+ aj Sv 12 ) (Sv 12 ) 

-2{ata 3 )(5v 12 )\ (7) 
the one-body density evolution (eq. (J5J) reduces to 

+ dB a ((afctjSvu) - (afaj) (Sv 12 )) 

+ dB* b ((5v 12 a+ aj ) - {aj aj ) (Sv 12 )) . (8) 
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It is interesting to notice that, although we consider a second-order perturbation theory for the residual interaction, 
the second order term exactly cancels out when approximation J7J) is used. Eq. (J7J corresponds to a gaussian 
approximation for quantal fluctuations. Therefore, eq. JHJ provides the stochastic equation of motion of one-body 
degrees of freedom associated to eq. (0 when neglecting part of the quantal fluctuations. The corresponding stochastic 
evolution of p reads 

dp=^r[h M F,p] + dB a {\ - p)U(p, a a )p 
ih 

+ dB* bP U'(p,a b )(l-p), (9) 

where Hmf denotes the matrix elements associated to the mean-field Hamiltonian while U(p,(T SL ) = 7 1 ? , 2(^ w i2( cr a)P2) 
and U'(p,a h ) = Tr 2 (p25vi2((Jb))- 

The stochastic one-body evolution given by eq. © contains also part of the information on correlations. Indeed, an 
approximate evolution of the two-body density, whose matrix elements are (af aj aia k ) = {kl | P12 1 ij) can be obtained 
through approximations similar to eq. J7J but preserving the symmetry of the two-body density: 



< afa+aia k vi2 > ~ < a+a+aia k X v 12 > + (< afa k v 12 X a+ai > - < a+a;i>i 2 >< a+a fc >) 
+ (< a+a fc >< a+a;wi2 > - < a+a; >< a+ a k v 12 >) 

- 2 (< afa k X a+ai > - < a+a ; >< a+a fe >) < v 12 > (10) 

and 

< a+ alj aia k v 12 vi2 > ~ < afa^a t a k X v\ 2 > +2 (< afa k Vi2 X a^a t vi 2 > - < afaivi 2 X a+a k vi 2 >) 

- 2 (< afa k X a+a; > - < a+a; >< a+a fc >) < v 12 > 2 . (11) 

Combining with Ito rules, the evolution of p±2 reduces to 

d < afafaidk >^ d{pupij ~ PkjPu) (12) 

indicating that the two-body evolution can be deduced from the stochastic evolution of p. Although eq. (|12fl is 
similar to the mean-field case, it contains correlations beyond mean-field. A similar situation occurs in the exact 
reformulation of self-interacting fermions with quantum jumps . 

In summary, the jump process described by eq. JSJ for both state vectors entering in D can be approximated by the 
jump process in one-body space given by eq. 10 if part of the quantal fluctuations are neglected. The advantage of 
this approximation is that expression JBJ for D is preserved along the stochastic path. In this work, we restrict ourself 
to this limit and eq. © will be referred to the incoherent Stochastic Mean-Field (SMF) dynamics. The properties of 
this diffusion process are described below. 

We consider that the single-particle states of |$ ) and initially verify 

((3 j \a i )=S ij . (13) 
Eq. @ can be simulated by quantum jumps for single-particle states given by 

|da 2 ) = + dB a (l - p)U(p, O] 

(14) 

Wi\ = (Pj\ [-§ h MF(p) + dB* b U'(p, <7b)(l - p)} . 

The latter quantum diffusion process has several attractive aspects. First, it can be easily verified that eq. 1131) is 
preserved along the stochastic path. Thus, the one-body density reads at all time p = \a.i) Consequently, the 
trace of the density is constant along the path: Tr(dD) = Tr(dp) = 0. In addition, p remains a projector, i.e. p 2 = p 
at all time. Finally, the total entropy S = —fcsTr(I) In D) is constant along the path. Indeed, since the density is 
given by eq. 10, S(D) identifies with the one-particle entropy S(p). Using equation and Ito rules, we obtain 
dS(p) = 0. Despite a constant entropy, the SMF induces correlations beyond mean- field. Indeed, starting from an 
initial two-body density p\2 — A{p\p2), after one time step, the average evolutions of the one- and two-body density 
matrices read 

dp = §[h M F,p] 

(15) 

d-plt = § [h MF {l) + W(2), pu] + dC 12 . 
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The labels "1" and "2" refer to the particle on which the operator is acting |3^|. dC\ 2 corresponds to correlations 
beyond mean-field associated to the stochastic one-body evolution given by eq. ©■ It reads 

dC 12 = ~^- {(1 - p x )(l - p 2 )U 1 {a a )U 2 {a a )p 12 

+Pi2U[(aM(a h )(l - Pl )(l - p 2 )} , (16) 

where the density dependence are omitted in U and U'. Eq. Ijl6(l clearly indicates that dC\ 2 is a second order term 
in perturbation. Note that it has a similar form as the second moment of the initial stochastic correlation used in ref. 
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FIG. 1: Top: Evolution of the root mean square radius (rms) as a function of time. Black circles correspond to the standard 
TDHF evolution while different lines correspond to different stochastic paths. Bottom: Error bars correspond to the rms 
evolution obtained by averaging over different paths while black circles correspond to the TDHF case. The stochastic simulation 
is performed for go = 500 MeV/fm. The average is taken over 200 trajectories. The width of the error bars correspond to the 
statistical fluctuations of the rms. 



B. Illustration of application 

To illustrate the SMF theory, we consider the monopolar vibration of a 40 Ca nucleus. The system is initially 
prepared in a pure state D = |$) ($|, where |<E>) is a Slater determinant solution of a constrained Hartree-Fock (CHF) 
equation. The CHF equation is solved assuming spherical symmetry and spin and isospin saturation. The Skyrme 
interaction of ref. [34[ is used in the mean-field. We assume in addition to the self-consistent mean-field, a monopolar 
constraint Xr 2 with A — 0.25 MeV.fm" 2 at t < fm/c [35|. At t — fm/c, the constraint is relaxed and two dynamical 
calculations are considered. The first one corresponds to the time-dependent Hartree-Fock (TDHF) evolution. In the 
second case, the SMF evolution described by eq. © is performed with a statistical ensemble of contact interactions 
defined by one stochastic variable, i.e. Sv\ 2 (a) = av res where v res is a contact interaction. In this case, U{a,r) takes 
the form U(cr,r) — agop(r), where go is a parameter measuring the strength of the perturbation. In both cases, 
evolutions are solved assuming spherical symmetry. 

The evolution of the root mean square radius (rms) obtained with TDHF is presented in figure ^ (filled circles) . 
The different lines displayed on the top part of figurenjcorrespond to the evolution of the rms along several stochastic 
paths obtained with go — 500 MeV/fm and a collision time t = 0.01 fm. In each case, the stochastic evolution 
differs significantly from the mean-field prediction. Bottom part of figure ^ shows a comparison between the TDHF 
evolution and the evolution of the rms obtained by averaging over the different stochastic trajectories. Interestingly 
enough, the average evolution identifies with the TDHF evolution. This example illustrates a special situation where 
the mean-field dynamics can be recovered from complex trajectories in many-body space |36| . However, significant 
fluctuations around the mean TDHF trajectories are observed. This is illustrated in figure [21 where the quantity 

A r = \J (r 2 ) 2 — (r 2 ) is displayed as a function of time for different values of go and A. We also computed as a reference 
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FIG. 2: Evolution of the dispersion of the rms as a function of time for different values of go- The different curves from bottom 
to top corresponds respectively to go = 100, 250 and 500 MeV.fm. Solid lines and open circles correspond respectively to an 
initial constraint A = MeV.fm -2 and A = 0.25 MeV.fm" 2 . 

the mean- field width omf-, corresponding to the quantal fluctuations of r 2 estimated in the mean- field approximation, 
leading to umf = 6.3/m 2 . As we do expect in the weak coupling approximation, the additional fluctuations A r 
induced by the stochastic term are much smaller than <jmf- It turns out that the dispersion is properly parametrized 
by the formula A r = Aq(1 — e~ Tat ) where Ao is proportional to go while Tq is independent of it. Therefore, while 
the average evolution of the rms collective variables is not affected by the stochastic process, fluctuations around the 
mean value increases and saturates as we do expect in Brownian motion. Interestingly enough, the behavior observed 
here is very similar to the description of a quantum oscillator 37] with Nelson stochastic mechanics [38| replacing 
(h/m) 1 / 2 by go- Assuming a single collective state and using similar techniques as in ref. [37|. analytical expression 
can be obtained for A r where Ao is indeed proportional to go while To depends only on the oscillator frequency. Note 
however that a complete understanding of the Brownian process presented here pass through the linearization of eq. 
@ as in ref. |l7j . 

C. Critical discussion 

In this section, it has been shown that the presented SMF theory can be applied to account for fluctuations beyond 
mean-field. The above example can only serve as an illustration due to the very schematic residual interaction used 
and to the very small time r (see estimation in ref. [3^. lioTp. Besides the simplicity of the force, it is important to 
note that the above theory only gives a partial answer to the simulation of correlations beyond mean-field for realistic 
nuclear systems. Indeed, as clearly seen in figure ffl. while fluctuations of one-body observables are increased, the 
average evolution of the rms identifies with the TDHF case. Therefore, no additional damping is observed in SMF in 
opposite to what is obtained in extended Time-Dependent Hartree-Fock (extended TDHF) [T3]- Indeed, dissipative 
aspects present in the memory kernel of extended TDHF are not included in the present framework. A careful 
analysis demonstrates that the absence of a collision term could be assigned to the approximations made on quantal 
fluctuations (equation Q, I|1U|) and 

Under these approximations, it is however possible to show that a fermionic system submitted to a statistical 
ensemble of residual interactions 5v 12 can be treated by a jump process related to its mean-field (equation ©)• The 
form of the noise is a second critical aspect for nuclear physics. In fact, it is expected that the residual interaction 
is dominated by the two-particle two-holes (2p-2h) channels. Due to the mean-field nature of the noise in eq. 
these components cancel out in the stochastic part. Therefore, while the above SMF can be of great interest, some 
important aspects for nuclear physics are missing. In the next section, we discuss the possibility to use quantum jump 
process in one-body space in a more general framework. 

III. STOCHASTIC MECHANICS WITH DISSIPATION 

In order to generalize the stochastic method described previously, we consider directly the evolution of the many- 
body density D associated to eq. Q. Using the short memory time approximation for 5v\2, the evolution of D can 
be recast as 

AD = ^-[Hmf, D] - ^[8vi2,[8vi2,D]], (17) 
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In the previous section, the term 5vi2DSvi 2 has been neglected due to the fact that the average evolution of wave- 
functions were directly considered. Due to this approximation, it was possible to have uncorrelated noise for |$ a ) and 
($f,|. Here, this contribution is not neglected. As a consequence, the uncorrelated noise assumption is not possible 
anymore. 

The one-body density matrix equation of motion associated to (|17|l is given by 

§ = ^[h MF (p),p}-^V(p), (18) 

where g is a real constant, g = r/h 2 . V(p), called "dissipator" hereafter, corresponds to the average effect of the 
residual interaction and reads 



(j\V\i) = {[[a+a 3 ,Sv 12 ],5v 12 ]). (19) 

We assume that the system is initially in a pure state described by a Slater determinant \$>(t)} formed of N orthonormal 
single particle states denoted by \a). The associated initial one-body density matrix reads p = ^ a \a) (a\. Having 
in mind the nuclear many-body problem, we assume that only 2p-2h components of Svi 2 are not equal to zero. 
Completing the hole states by a set of particle states, denoted by \ce), we have: 

Sv 12 (t) = 44 ("P <*p) OaO/3- (20) 



\ 
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V can then be recast as 



where F\ 2 is equal to 



V(p)=Tr 2 [v 12 ,F 12 ], (21) 



F12 = - {(1 - pi) (1 - p%) V 12 pip 2 

- pip 2 vi2 (1 - pi) (1 - p 2 )} ■ (22) 

Expression (|21|l takes a form similar to the collision term generally obtained in extended TDHF The dissipator 
V(p) can be further transformed. Indeed, 5v\ 2 given by eq. lf2*U|l can always be decomposed as (see for instance |41f ) 

n 

where X n are real and the O n correspond to a set of commuting Hcrmitian one-body operators written as O n = 
Efia (a \O n \ a) a^a a . Reporting in eq. ff^Tf) . T>(p) can be recast as 

V(p) = T m „ [O n O mP + P O n O m - 20 m pO n ] . (24) 

mn 

The coefficient T mn are given by 

r m „ = -A m A„Tr(O m (l - p)O n p). (25) 

We recognize in this expression, the quantum covariance between the operator O n and O m , i.e. Tr(O m (\ — p)O n p) = 
(O m O n ) — (O m ) (O n )- Expression l|24|l has the form of the dissipator appearing usually in the Lindblad eauation[25|. 
Therefore, we have shown that the evolution of one-body degrees of freedom associated to equation (11711 identifies 
with a Markovian quantum master equation generally obtained in quantum open systems. A large amount of work 
is devoted to the simulation of such master equation by quantum jump methods (see for instance |25l l42l I43L I44I |45| ) 
and one can take advantage of the most recent advances in this field. This aspect has however rarely been discussed in 
the context of self-interacting system. In the following, the associated diffusion process is precised and we show that 
it indeed corresponds to jumps between Slater determinants. The stochastic Schroedinger equation for single-particle 
wave-function is finally given. 
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A. Explicit form of the stochastic process 

Following ref. pEj . we introduce the Hermitian positive matrix T with components T mn . An economical method 
to introduce quantum jump process [25| is to use the unitary transformation u that diagonalizes V, i.e. Y — u^ 1 ju, 
where 7 is the diagonal matrix of the eigenvalues of T. New operators A k can be defined by the transformation 
A k — J2 n u knOn- The dissipator is then recast as 

V(p) = [Alp + pA\ - 2A kP A k ] . (26) 

k 

The last equation can be simulated using the average over the stochastic mean-field dynamics: 

dp = % [h MF (p),p] ~ g%V{p) + db sto , (27) 
in 2 

where db B to is a stochastic one-body operator which, using Ito rules |32j |. reads 

db sto = 5^{dWk(l-p)i4*p + dWJpA fc (l-p)}. (28) 

k 

Here dW k denotes stochastic variables given by dW k — ~id^ k ^/g^ k: where d£ k corresponds to a set of real gaussian 
stochastic variables with mean zero and dt; k d£ k i = 5 kk idt. 



B. Nature of the stochastic process in Hilbert space 



It is worth noticing that the proposed dissipative equation and its stochastic counterpart are only well defined if the 
density is initially prepared as a pure Slater-determinant state. Wc now turn to the essential properties of equation 
(l27jl . First, it preserves the number of particles Tr(dp) = 0. In addition, if initially p 2 = p, then 

dpdp - gj [pV(p) + V(p)p] = -gjV(p) (29) 

which is obtained using Ito stochastic rules and retaining only terms linear in dt. The last expression demonstrates 
that {p + dp) 2 = p + dp. Thus, p remains a projector along the stochastic path. As a consequence, the pure state 
nature of the many-body density matrix is preserved along the stochastic path, i.e. D — |$(t)) (&(t)\ where |$) is a 
normalized Slater determinant at all time. The associated stochastic Schroedinger equation for single-particle states 
reads 



d\a) = \^-h MF {p) + ^dW k {l- p)A k -g-Y,lk [A\p + P A kP A k - 2A kP A k ] j \a) . (30) 

This last expression can be directly used for practical applications. 

In this section, assuming that the residual interaction can be written in terms of 2p — 2h components, we have 
shown that the dissipative dynamics of one-body degrees of freedom can be simulated by quantum jumps. The 
stochastic method differs significantly from the simplified scenario considered in previous section. First, the above 
theory does not require to introduce generalized many-body density matrix formed of two Slater determinants since 
D = ($(t)l along each path. As a counterpart, the numerical effort required to treat the new SMF which 

includes dissipative aspects is significantly increased. Indeed, in section Til Bl we have shown that fluctuations can be 
introduced using simplified residual interaction and a limited number of noise terms. In the new SMF, the numerical 
effort is directly proportional to the number of one-body operators entering in eq. I|23|l . This situation is similar to 
stochastic methods used in nuclear structure studies |19| . In that case, large numerical effort is required. For instance, 
if we assume that a physical system is described in a mesh, then the number of noise terms is is a priori as large 
as the number of mesh points. The numerical implementation of the quantum jump process in its full complexity is 
expected to be still difficult with present computers capacities. Therefore, specific numerical techniques as well as 
truncation procedures should be developed in order to implement the method in many-body dissipative dynamics. 
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IV. CONCLUSION 

In this work, we have presented a discussion on the possibility to replace the dynamics of interacting fermions by 
quantum jumps. Assuming a weak coupling approximation and a short memory time for the residual interaction, 
two scenarios have been considered. Focusing on one-body degrees of freedom, two different approximations lead 
to equations of motions for the one-body density matrix that can be treated by quantum jumps between Slater 
determinants. 

In the first stochastic mean-field process, part of the quantal correlations between p and 5v\2 are neglected along 
the path. The SMF is illustrated in the monopolar vibration of a calcium nucleus. In this case, while expectation 
values of one-body observables are unchanged, fluctuations are increased compared to mean-field. In the presented 
application, the residual interaction is rather simple. However, more complex statistical ensembles can be used like 
two-body random interactions. This method suffers from the absence of dissipative effects and appears to be too 
simplified for the nuclear many-body problem. 

In the second part of this work, we show that the dynamics of correlated systems can be described by a quantum 
jump process even if no approximation on the quantal correlations are made. Restricting to 2p-2h residual interactions, 
a guideline is given for transforming the evolution of the one-body density matrix evolution into a Lindblad equation. 
The stochastic process corresponds to quantum jumps between pure state densities D — |$) ($| where |4>) is a Slater 
determinant. The associated stochastic equation of motion for single-particle wave-functions is given. 

Finally, we would like to mention that the presented framework does not account for non-Markovian effects. It has 
however been shown in ref. 15] that the memory effect might be important in the nuclear context. Promising work 
are actually devoted to incorporate non-Markovian effects in quantum Monte-Carlo methods |46l EtIEsI Efij . 
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